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We study critical point finite-size effects on the behavior of susceptibility of a film placed in 
the Earth's gravitational field. The fluid-fluid and substrate-fluid interactions are characterized by 
van der Waals-type power law tails, and the boundary conditions are consistent with bounding 
surfaces that strongly prefer the liquid phase of the system. Specific predictions are made with 
respect to the behavior of ^He and ^He films in the vicinity of their respective liquid-gas critical 
points. We find that for all film thicknesses of current experimental interest the combination of van 
der Waals interactions and gravity leads to substantial deviations from the behavior predicted by 
models in which all interatomic forces are very short ranged and gravity is absent. In the case of 
a completely short-ranged system exact mean-field analytical expressions are derived, within the 
continuum approach, for the behavior of both the local and the total susceptibilities. 
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I. INTRODUCTION 

In the case of the liquid-gas critical point, gravitational 
effects become important as the isothermal compressibil- 
ity (i.e. the susceptibility) increases to a large value in 
the vicinity of that point and, indeed, diverges as that 
point is approached. This leads to a vertical gravity- 
induced density gradient that grows as the critical point 
is approached. Even in a relatively small experimental 
cell with a vertical dimension of 0.5 mm, filled with ^He 
at its critical density, the density stratification between 
the cell top and bottom is approximately 6% at a re- 
duced temperature of t = 10~^ [T]. Only the density at 
the middle of the cell remains at its critical value. In 
order to interpret precision critical point measurements, 
one must develop models that incorporate the effects of 
gravity. A number of theoretical studies have been per- 
formed on the effect of gravity on measurements near the 
superfluid transition in ^He. Models of the ^He specific 
heat in bulk and finite-size samples near the superfiuid 
transition have been successfully tested by high precision 
measurements in ground-based laboratories jSJ [3] and in 
microgravity [4| [5] . 

In a previous article [6 , the authors of this paper re- 
ported the results of a theoretical calculation of the ex- 
pected finite-size effects on the isothermal susceptibility 
near a liquid-gas critical point in a zero gravity envi- 
ronment. That investigation was performed for a thin 



film between surfaces that both strongly prefer the liq- 
uid phase. It is hoped that predictions of this study will 
be tested by future experimental studies in space [7 . In 
the present paper, we extend those calculations including 
the effects of gravity. The ^He and ^He liquid-gas crit- 
ical points were chosen for this theoretical investigation 
because these systems are devoid of impurities and many 
bulk thermophysical properties have been measured near 
their liquid-gas critical point. 

In this article we will discuss the behavior of the sus- 
ceptibility of a film of a non-polar fiuid of, say, ^He 
or ^He, having a thickness L in which the intrinsic in- 
teraction is of the van der Waals type, decaying 
with distance r between the molecules of the fiuid as 
^-(d+cr) Yiere d is the dimensionality of the system 
while cr > 2 is a parameter characterizing the decay of 
the interaction. The film is bounded by a substrate, say 
Au plates, that interacts with the fiuid with similar van 
der Waals type forces, i.e. of the type J^'^ ~ z~^% where 
z is the distance from the boundary of the system while 
cTg > 2 characterizes the decay of the fluid-substrate po- 
tential. For realistic fluids d = a = cFs =3. 

According to flnite-size scaling theory P| [8| l9l , the be- 
havior of the susceptibility in a fllm of a fluid placed 
in an external gravitational fleld, governed by dispersion 
forces and subject to (+, +) boundary conditions, i.e., 
conditions that strongly favor the liquid phase of the fluid 
over the gas one, is 



X{t, All, L) - Xbulk(t, Am) - Kface(*> Am) + X^uXe A/.)] = (1.1) 

L^'-X{L/^u a^(/?AM)i^/^ ag{pg)L^+''/\ h^,,L-'^%hL-"\a^L-'^), 

I 

where Xbuik(i, A/i) is the bulk susceptibihty, Xsurface(*' A^) and XsurfacT(*' Am) are the sur- 
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face susceptibilities — the result of the exis- 
tence of two surfaces bounding the system — 
6(T) = ioo{T ^ T+,A/i = 0) ^ i^\t\-^ is the 
bulk correlation length, t = (T — Tc)/Tc is the re- 
duced temperature, Tc is the bulk critical temperature, 
A/i = /i — /ic is the excess chemical potential, while 
lie is the bulk critical chemical potential, g is the 
external gravitational field, ^q", a^^ag^ h^^s^ b and are 
nonuniversal metric factors, and 

uJs = (Ts-{d^2-r])/2, uJb = a^r]-2. (1.2) 

The quantities i/. A, and uo are the universal critical 
exponents for the corresponding short-range system. 

With respect to their bulk critical behavior, the nonpo- 
lar classical fluids belong to the so-called Ising, or 0(1), 
universality class. When d = 3 this universality class is 
characterized by critical exponents \iO^ 



T] = 0.034, 7 = 1.2385, u = 0.631, 



and 



a = 0.103, = 0.329, 6> 



0.53. 



(1.3) 



(1.4) 



As we have already stressed, the methods and ideas to 
determine the behavior of the susceptibility if thin films 
of non-polar fluids used in the current article are general 
and can, in principle, be applied to any non-polar fluid; 
however, but we will exemplify these methods via the 
study of ^He and ^He films. 

A central question pertains to the extent to which the 
critical exponents listed above describe the experimen- 
tally observed behavior of these substances near their re- 
spective liquid-vapour critical points. In the early years 
of the development of critical behavior, there were at- 
tempts to experimentally measure the critical exponents 
of the He isotopes to be compared with theoretically pre- 
dicted values; see, e.g., [IT]. However, it became clear 
that the asymptotic power-law region is very limited in 
ground-based measurements and accurate analyses must 
take into account correction-to-scaling terms as well as 
gravity effects; see, e.g., [12 . These issues are partic- 
ularly important in the case of the He isotopes, which 
have the largest gravity effect (see, e.g., the discussion of 
Table I in [7]). Because of these problems, most recent 
experimental measurements, particularly those in the He 
isotopes, have been compared to models using the the- 
oretical critical parameters obtained from those models; 
see[13 . In light of this history, we will utilize theoreti- 
cally derived values of the relevant critical exponents. 



From (1.2) and with a = (jg =3, as in physical van 



der Waals interactions, one has cjg = 0.52 and cj^ = 1.03. 
Since uoi > uo > uOs > for L large enough., one can 
expand the r.h.s. of Eq. with the result that the 

leading finite-size behavior of the susceptibility near the 
bulk critical point is given by the properties of the cor- 
responding short-ranged system. This is true because 
all the dependences on the long-ranged tails of the van 



der Waals interaction in are reflected via the fac- 

tors 6, hw^s and 6 is proportional to the strength 
of the fluid-fluid interaction [IKTS], hy^^s reflects the 
contrast between the fluid-fluid and the substrate-fluid 
J^'* effective interaction at Tc [6l[8] (see below), while the 
field a^y, associated with the Wegner-type corrections to 
scaling, in general incorporates contributions due to the 
long-ranged tails of the interaction [16 . In [6l [8] it is 
demonstrated that such an expansion is admissible only 
when L is much larger than some critical thickness Lcrit, 
i.e., when 



L » ierit = {2-+^\h^Ay'^ - 200 Co+ \h 



11.918 



A. 

(1.5) 

For most systems ^q" is of the order of 3 A and h^^s is 
of order of 1. For some systems, e.g., like ^He and ^He 
bounded by Au the dimensionless constant h^^s can be 
as large as 4 [6l|8]. The constraint (1.5) represents the 
^'relevance-irrelevance^^ criterion for the van der Waals 
forces with respect to the behavior of finite-size quanti- 
ties in van der Waals thin films; when L ^ Lcrit such 
forces can be neglected, while when L < Lcrit they must 
be taken into account in, say, the determination of the 
behavior of the finite-size susceptibility. One can also 
formulate a criterion for the relevance of gravity. From 
Eq. ( |1.1| ) it is clear, that gravity is a relevant variable. 
If, however, {pg)L^^^~^^ <C 1 the influence of gravity can 
be neglected and will play no essential role in the be- 
havior of any finite-size quantity. In the opposite case 
its role is crucial and must be taken into account. Note 
that this criterion also connects the relevance of gravity 
to the thickness L of the films; for thin films it is negligi- 
ble, while in the case of sufficiently thick films it is not. 
In the case of ^He and ^He we will discover that films 
with L = 1000, 2000 or 4000 liquid layers are, in this 
sense, thin films, while a film with, say, L = 8000 layers, 
can be considered thick for the purposes of assessing the 
influence of gravity on critical point behavior. 

In this article we discuss the finite-size behavior of 
the susceptibility of van der Waals fluid films bounded 
by two flat substrate plates situated perpendicular to 
the Earth's gravitational field (i.e. horizontal), both of 
which strongly prefer the liquid phase of the system. The 
schematic phase diagram of such a system in the (T, 
A/i) plane is shown in Fig. [l] The solid line A/i = 0, 
T < Tc represents the bulk gas-liquid phase coexistence 
line. The liquid-gas coexistence curves A/icap correspond 
to capillary condensation transitions for L = Li, L = I/2, 
and L = Ls where Li < L2 < T3. The case ^ = 
has been extensively studied, and the picture presented 
here is in accord with the results of Refs.[T7l [HI [T9l [20]. 
When no gravity is present within the system, the fact 
that A/ic,Li < simply expresses the preferences of the 
identical walls (see the cases Li and L2). In the pres- 
ence of a gravitational field orthogonal to the bounding 
plates, the action of gravity pulls the molecules of the 
fluid away from the upper plate, leading to a more gas- 
like phase near that plate and a more condensed liquid- 
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FIG. 1: Color online. The schematic phase diagram of a d- 
dimensional film system for various thicknesses L in the pres- 
ence of a gravity subject to boundary conditions that strongly 
favor the liquid phase at the plates bounding the fluid. 



nor van der Waals type interactions are present, a per- 
fect data collapse can be achieved, see Fig. [3) where the 
data are plotted in the same way as in the current fig- 
ure. Note also that both the van der Waals interactions 
and the presence of gravity reduces the magnitude of the 
susceptibility. This is due to the ordering effect of the 
van der Waals interactions and gravitational field. In the 
cases L = 1000, 2000 and 4000, the van der Waals in- 
teractions constitute the important influence leading to 
the effect, and this reduction is approximately given by 
a factor of 2. On the other hand, when L = 8000 the 
reduction in the maximum value of the susceptibility is 
by a factor of 3 times and is primarely due to the influ- 
ence of gravity. We also note that the maximum of the 
susceptibility as a function of Aja also changes its loca- 
tion; while for L = 1000, 2000 and 4000 the maximum 
occurs when A/i < with x)^^^ = (/3A/i)L^/^ = 0(1), 
for L = 6000 and 8000 the maximum is at A/i > with 
x^^^ ^ /3cgL/2 ^ 1. One of the aims of the current 



like phase region near the lower one. The larger L and 
^, the stronger this effect. Thus, for L large enough (see 
the case with L = Ls) the critical point of the finite sys- 
tem lies above the A/i = line, i.e. at Afic,L > 0, which 
stabilizes the liquid phase of the fluid. Away from the 
critical region, the shift in the phase boundary relative 
to the bulk coexistence line A/i = is proportional to 
while within the critical region it is proportional 
to where A and v are the standard bulk critical 

exponents. The lines of first-order phase transitions end 
at {d — l)-dimensional critical points Tc{Li) with coor- 
dinates {Tc^L-, Afic,Li)^ ^ = 1,2,3, the positions of which 
vary with L and depend on the presence of gravity ^, 
as well as on the presence and on the strengths of the 
fiuid-fluid and the substrate-fluid interactions. 

For large L these points are located close to the bulk 
critical point Tc with coordinates {Tc^ Aji = 0): Tc^l — 
Tc - L-i/^ and A/ic,L - ^f^gL - L"^/^. Since the 
fluctuations in systems of reduced size are stronger, one 
typically has Tc,l^ < Tc. In part, the structure of this 
phase diagram is reflected in Fig. [2j where the behavior 
of the finite-size susceptibility in thin films with thick- 
ness L = 1000,2000,4000,6000 and 8000 layers is shown 
as a function of the scaling variable = (/3A/i)L^/^ at 
the bulk critical temperature T = Tc of the correspond- 
ing infinite system (with ^ = 0). Both the presence of 
the van der Waals interaction between the fluid parti- 
cles and between the substrate and the fluid, as well as 
the gravitational field of the Earth are taken into ac- 
count. One observes a clear lack of data collapse. For 
relatively thin films— with L = 1000, 2000 and 4000— 
this is due to the role of the van der Waals interaction, 
while for relatively thick films with, say, L = 8000 the 
absence of data collapse is due to the presence of gravity. 
The L = 6000 case illustrates the intermediate situation 
when the influence of van der Waals interactions fades 
away and the gravity steps in as a factor mainly respon- 
sible for the lack of data collapse. When neither gravity 
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FIG. 2: Color online. The behavior of the finite-size suscepti- 
bility X in thin films with thickness L = 1000, 2000, 4000, 6000 
and 8000 layers at the bulk critical temperature T = Tc of the 
corresponding infinite system as a function of the scaling vari- 
able — {f3Afi)L^^^ . The corresponding scaling variable 
that governs the dependence on the gravity is proportional to 
Xg ^ {l3g)L^^^^^ and, thus, the gravitational effects gradu- 
ally set in with increase of L in the behavior of the finite-size 
susceptibility. 

article is to explain how the above curves have been ob- 
tained and to elucidate why, as we believe, these curves 
ought to resemble the ones obtained in experiments with 
real liquid-gas systems. In order to facilitate contact with 
experiment we will, in addition to presenting the behav- 
ior of X ill the (T, A/i) plane, derive the corresponding 
dependence of x iii the (T, Ap) plane, where Ap = p — pc, 
with pc the critical density of the fluid. 

The structure of the article is as follows. First, in Sec- 
tion [n] we present a precise formulation of the model of 
interest. The corresponding simplification of the model 
in the case of a film geometry and the analytical expres- 
sions needed for its numerical treatment are presented in 
Section [nil The results for the behavior of the finite-size 
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FIG. 3: Color online. The behavior of the finite-size suscepti- 
bility X at T = Tc in thin films with short-ranged interaction 
and with L = 1000,2000,4000 and 8000 layers. Note that 
the maximum of the susceptibility is at negative values of 
— {PAiii)L^^^ Note the excellent scaling in the behavior 
of the curves for the different L's, contrary to the behavior of 
the same system but with gravity and van der Waals interac- 
tions with the bounding the system plates taken into account 
- see the previous figure. 



susceptibility at the critical point as a function of /i, L 
and g are presented in Section [IV| while the correspond- 
ing results for T > Tc and T < Tc are given in Section 
[V| The article closes with a discussion and concluding 
remarks. 



II. THE MODEL 

Following Refs. [6] and [8] we consider a lattice-gas 
model of a fluid confined between two parallel flat plates 
at a distance L with a grand canonical functional r^[p(r)] 
given by 



mr)] = kBTj2{p{r)\n[p{r)] 
rec 

+[l-p(r)]ln[l-p(r)]} 

+ ^[FO'^)(z)+^z-/i]p(r). (2.1) 
This expression is to be minimized with respect to the 



local number density p(r). The functional (2.1) is the 
simplest model that captures the basic features of sys- 
tems with both van der Waals interactions and gravity 
taken into account. It can be viewed as a modification 
of the model utilized by Fisher and Nakanishi in their 
mean-field investigation of short-range systems [171 [IH] 
and in th e ab sence of gravity. 

InEq. (2.1) w^^\r-r^) = -47*^^) (r-r') is the non-local 



a simple cubic lattice in the region occupied by the fluid. 
We consider a fluid system with geometry oo^~^ x [0, 
where the region < 2: < L is occupied by fluid. Here 
and in the remainder of this paper, all length scales are 
taken in units of the lattice constant a, which is of the or- 
der of a molecular diameter — which means that length is 
expressed as a dimensionless quantity — so that the par- 
ticle density p(r) is dimensionless and varies within the 
range [0, 1]. In Eq. (2.1), the terms in curly brackets cor- 
respond to the entropic contributions, while the other 
terms are directly related to the interactions present in 
the system. The term proportional to g reflects the pres- 
ence of gravity. It is assumed that the gravitational field 
is along the z direction, i.e. is perpendicular to the planes 
bounding the fluid. The external potential V^^'^\z) re- 
flects the interaction between the molecules of the fluid 
and of the constituents of the two substrates bounding it. 
For an individual wall V^^'^\z 00) ~ z~^ with a = 3 
for a genuine van der Waals interaction. In the current 
treatment we will assume that the bounding substrates 
are identical on the both sides of the film and that they 
strongly prefer the liquid phase of the fluid. In Eq. (2.1), 
jii is the chemical potenti al. 

The variation of Eq. (|2.1[) with respect to p{r) leads 



to the equation of state for the equilibrium density p* (r) 



2p*{r) - 1 



tanh 



^w;(')(r-r')/9*(r') 



(2.2) 



The advantage of this equation is that it lends itself to 
numerical solution by iterative procedures. For a partic- 
ular geometry and surface potential V'^^'^^z) the solution 
determines the equilibrium order-parameter profile p* (r) 
in the system. Inserting this profile into Eq. (2.1) one 



obtains the system's grand potential. To avoid the dou- 
ble sum in Eq. (2.1 ), which is inconvenient in a numerical 



treatment, we make use of the relationship below, which 
is easily derived from Eq. (2.2) 



I ^ w{r-v')p*iv)p*{r') 

]^Y[^,-Viz)-gz]p*{r) 
kBTYp*{r) arctanh[2/j*(r) - 1], 



(2.3) 



_ 1 
2 



which, when inserted in Eq. (2.1), yields 



rec 



kBT{p* {r)\n [p*{r) 



coupling between the constituents of the fiuid, while C is 



-[1 - p*(r)] In [1 - p*(r)] - p*(r) arctanh[2p*(r) - 1]} 

1 
2 



[p-V{z)-gz] p*(r) 



(2.4) 
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Note that in ( |2.4| ) p*(r) is no longer a free functional 
variable, but is the solution of Eq. (2.2). 

Denoting 0*(r) = 2p*(r) — 1 and A/i = /i — /i^ , where 
jjLc = \ X^r' '^(^ ~ ^')^ equation of state (2.2) can be 
rewritten in the standard form 



0*(r) = tanh 



/3^J(r,r')r(r') 



(A/x - AT/(z) - gz) 



(2.5) 



where J(r — r') = —w(y — The bulk properties of 

the model are well known (see, e.g.,[22l|23] and references 
therein). We recall that the order parameter ^* of the 
system has a critical value (/>* = which corresponds to 
Pc = 1/2 so that 0* = 2(p* - pc). The bulk critical point 
of the model is given by {13 = Pc = [Zlr "^(^)]~^' ~ 
jic = ~2 J]]^ J{^)) with the sum running over the whole 
lattice. Within the mean-field approximation the critical 
exponents for the order parameter and the compressibil- 
ity are jS = 1/2 and 7 = 1, respectively. As has been 
shown in [6l[8], the surface potential AV(z) is 

AF(z) = 5vs [{z + 1)-" + (L + 1 - z)-"] , (2.6) 

where 1 < z < L— 1, and where contributions of the order 
^, etc. have been neglected, the quantity 



6v. 



_4^((i-l)/2 



{psf^'-Pcf) (2.7) 



ar(^) 

is a (T- and /i-independent) constant, the quantity 

J(r) = J7(1 + |r|'^+-), (2.8) 

is a proper lattice version of —w{y)/A as the interaction 
energy between the fiuid particles, and 



J^'^(r) = J^'^ 



\d+(T 



(2.9) 



is the interaction between a fiuid particle and a substrate 
particle. Here ps is the number density of the substrate 
particles in units of a~^. Note that the effective potential 
5vs is the result of the difference between the relative 



strength of the substrate-fiuid interaction for a substrate 
with density ps and that of the fiuid -fiui d interaction for 
a fiuid with a density p^ In Eq. (2.6), the restriction 



z >1 holds because we consider the layers closest to the 
substrate to be completely occupied by the liquid phase 
of the fiuid, which implies that we consider the strong 
adsorption limit, i.e., /?(0) = p{L) = 1; therefore the 
actual values of AV{0) = AV{L) will play no r ole. 
In terms of the quantity ^ the functional (2.1 ) becomes 



1 - cPir) 



1 + ,^(r) 



In 



(r) 



-lY.[Aii-AV{z)-9z]<l>{r) (2.10) 



rGC 



I J2 J(r,r')(/.(r)(/.(r') + aeg, 



r,r'e£ 



where 



nr. 



-E 



A/i - AV{z) - gz ■ 



E 

r'eC 



J(r,r') 



(2.11) 

does not depend on <p and therefore is a regular back- 
ground term. An expression similar to the one in Eq. 
(2.4), which avoids the double sum and thus is more con- 



venient for numerical procedures, can also be obtained. 
With the identifications (j){r) ^ m(r) and 

H{z) = ^[Afi-AV{z)-gz] (2.12) 

one can rewrite the above expression for l][^(r)] as a 
functional of the effective magnetic density m(r): 



F[m(r)] = {n- n,eg), 



(2.13) 



which defines the free energy of a magnetic system at 
temperature T and in the presence of an external local 
and spatially varying magnetic field H{z). Explicitly, one 
has 



/Qz^r ^ \^ / ^ + i [1 + m(r)] 1 - m(r) ri-m(r)]\ v^./n / n ^ t^/ /a a a a /^ 
/3F[m(r)] = 2^ I In + In | - 2^ /i(z)m(r) - - 2^ i^(r, r )m(r)m(r ), 

I ^ 

where i^(r, r') = /3J^(r, r') is the non-local coupling be- In the remainder of this article we shall make use of 

tween magnetic degrees of freedom, h{z) = /3H{z) is an this connection between the fiuid and magnetic systems 

external magnetic field and the magnetization m(r) is to in order to exploit existing theoretical results for both of 

be treated as a variational parameter. them. 
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III. FINITE-SIZE BEHAVIOR OF THE MODEL 
IN A FILM GEOMETRY 



In [6] it has been shown that the function J{z) can be 
written in the form 



We win be interested in a system with a film ge- 
ometry. Because of the symmetry of the system one 
has 0(r) = 0(r||,z) = (j){z)^ where r = {r||,z}, i.e., 
the order parameter profile {^(z),0 < z < L}^ with 
0(0) = 4>{L) = 1, depends only on the coordinate perpen- 
dicular to the plates bounding the van der Waals system. 
In this case Eq. (2.5) becomes 



0*(z) = tanh 
where 



(3j2j{z-z')cP*{z')+h{z) 



(3.1) 



J{z) ^ -^(^11 - , ^) = E -^(^11 ' (3.2) 



J{z) = f[c^_,S{z) + cZi[S{z-l) + S{z + l)] 

+gd{z)e{z - 2)] , (3.3) 



where 5{z) is the discrete delta function, while 0{z) is the 
Heaviside function. Explicitly, for d = a = 3 one has [6] 



5 1 + Inl' 



3.602, 



(3.4) 



C"" = TT 

^ 3 



(-l)i/3^o( V 2 - 2iV37T) - (-l)2/3Xo( V 2 + 2iV37r) 



3 V\/3 



and 



^3(2;) 



\/3 arctan ^ ^ ^ —ln(l+ +^ln(l 



ln2 « 1.183, (3.5) 



(3.6) 



where Ko{x) is a modified Bessel function of the second 
kind. Following [6] and taking into account the presence 
of gravity in the system we study, the layer magnetic field 
h{z) is 



h{z) 



where 



1 



{z^iy 



zY 



(3.7) 



1< z < L - 1, 



the model under consideration. In order to have within 
the current mean-field model the same order of the finite- 
size effects due to hy^^s as in real systems we take in our 
model calculations cfq = 3.5. This value will be used in 
the remainder of the article whenever the substrate-fiuid 
interaction is taken into account. 

With respect to the behavior of the total susceptibility 
X of the system per unit particle it was shown in [6 that 



(3.9) 



(3.8) 



refiects the relative strength of the fiuid-wall and fluid- 
fiuid interactions, respectively. The above expression 
takes into account the fact that the substrate occupies 
the region R^-^ x [L + 1, 00] U R^"^ x [-{L + 1), -00]. 
For ^He and ^He bounded by Au surfaces in it has 
been shown that h^s = 4. Note that (3.7) is derived for a 
system with a = ag. According to finite-size scaling the- 
ory the finite-size effects due to the surface field h^^s are 
controlled by For = a = a, =3, an 

Ising-like system, this leads to h^^g/VL^ where the value 
of 77 = 0.034 has been neglected, i.e., 7^ = was used. Let 
us recall that within a mean-field treatment with respect 
to the critical behavior, the effective spatial dimension 
is (i = 4 irrespective of the actual spatial dimension of 



where R is the inverse matrix of the matrix R with 
elements 



■pj{z-z'), 



(3.10) 



while the "local" susceptibility, which reflects the re- 
sponse of the system from a given layer is 



xK^) = E(i^-') 



(3.11) 



with 

Xi{z) = EG(^,z*) (3.12) 

Z* 

= E(5(0,z)5(r|,^*)) - (5(0,z))(5(r|,z*)). 
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Obviously, x = E.XiW/(i + l)- 

For a fluid confined to a film geometry the natural 
quantity to consider is the excess grand potential nor- 
malized pe r uni t area A: Auj = \im.A^oo{^ — f^reg)/^- 
From Eq. (2.10) and (3.1) one obtains 



1 + <Piz) 



In 



1 + <Piz) 



iz) 



In 



1 - <j>iz) 



1 ^ 

-J2[^l^-^V{z)-gz]<p{z) (3.13) 

^ z=0 



z=0 z'=0 



The order parameter profile of the system is the one that 
provides the minimum of the above functional. For such 
a pr ofile {^*(z),0 < z < L}, which is a solution of Eq. 
(3.1), the above expression can be simplified to 




(3.14) 



--h{z)r{z) - - r{z) arctanh (^)] |, 

which is much more convenient for numerical evaluation 
since it does not involve the double summation present 
in dslsl). 



Eqs. (3.1), (3.2)-(3.10), and (3.14) provide the basis 



the susceptibility of a system in which both the van der 
Waals interaction between the molecules of the fluid and 
between the fluid and the constituents of the substrate, 
as well as the Earth gravity are taken into account. The 
procedure is as follows. First, we determine the order pa- 
rameter profile {^*(z),0 < z < L} hy solving iteratively, 
using the Newton-Kantorovich method, Eq. (3.1). How- 
ever, the solution of this equation depends, for a given 
range of parameters T and A/i, on the choice of the initial 
state of the order parameter profile. The two basic ini- 
tial states of the profile are i) di liquid-like state in which 
all the sites of the lattice are occupied by a particle, i.e. 
the state {(l){z) = 1,0 < z < L}, and ii) a gas-like state 
{0(0) = 1,0(L) = l,0(z) = 0,1 < z < L-1}. Thus 
one needs to calculate the profile starting from both of 
the two initial states. If the two final states coin cide t hey 
provide the unique minimum of the functional (3.13). If 



they differ one has to check which one provides the ab- 
solute minimum of the grad canonical potential. The 
simplest way to clarify that question is to calculate /3Au; 
viaEq. Kl^. 



for our numerical treatment of the finite-size behavior of 



The standard Ginzburg-Landau equation follows from 
Eq. (3.1), for small (/>, after taking into account that 
arctanh (0) + ^^/3 + 0(0^). A continuum version 
of the equation follows from the replacement 0(z — 1) + 
(p{z + 1) ^ 20(z) + (i)^^[z]. Obviously such a continuum 
version can also be constructed for the long-range system 
by adding the terms contributed by the function G{x), 
which is, in this case, a continuous function. Note that 
the function Q{x) is well defined everywhere for x > 
and not only for x > 1 as we actually need it in the 
lattice formulation of the theory. Thus, in the continuum 
formulation of the theory the integration can be extended 
over the region z G [0, L]. This does not change the long- 
range behavior of the magnetization profiles. Thus, in 
the continuum case the equation for the order parameter 
profile reads 



riz] + ^{r[z]r = h[z]+Ki c^riz) + cr 



2<P*{z)- 



if<j)*{z) 
dz"^ 



+ / g{\z-z'\')r{z')dz'}, (3.15) 



where K = jSJ'- and h{z), for < ^ < i, is 



The model with purely short-range interactions for 

A/i = and g = Q 



In the case in which all the interactions in the system 
are s hort-r anged and in the absence of gravity the equa- 
tion (3.15) for the order parameter profile in a continuum 
X 1 system can be written in the standard form 

h{z) = -/3Afi - -/3gz + hu,,s [z-^^ + {L - z)""^'] . 

(3.16) 

In continuum theory one defines the (+, +) boundary 

conditions via 4>*{0) = 4>*{L) = oo. - -j^ + a4> + u<p = h, (3.17) 
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where 



K 



and h 



h 



with 



1/{C2+2CT). 



(3.18) 



(3.19) 



The equations for the layer response function Xi{^)^ 
where 



dh 



h=0 



then reads 



+ (d + 3u4>^)xi = 1. 



(3.20) 



(3.21) 



Because conditions are identical at both bounding sur- 
faces of the system, the solutions of the above equations 
have to satisfy (t)'{L/2) = and x'{L/2) = 0. Under the 
(+, +) boundary conditions envisaged here (strong ad- 
sorption) one has, in addition, 0(0) = 0(L) = +oo, and 
X(0) = x{L) = 0. 

When h = the magnetization profile is known exactly 



a) when Xt = aL^ > — tt' 



<P{z) = L 



2K{k) 



dn[2K{k)C ; k] 
sn[2K{k)C 



2K{k)C 



dn[2i^(fc)C; k] 



, (3.22) 



sn[2K{k)C;k] 
where fc^ > is to be determined from 

Xt = [2K{k)f{2k'^ - 1), (3.23) 
and C = z/L, i.e., < C < 1- 
b) when Xt < — tt^ 



V u \sn[2K{k)C;k] 



where > is to be determined from 

Xt = -[2K{k)]\P^l), (3.25) 

and C = z/L, i.e., < C < 1- 

Here K{k) is the complete elliptic integral of the first 
kind, dn(C; k) and sn(C; k) are the Jacobian delta ampli- 
tude and the sine amplitude functions, respectively. The 
bulk critical point T = Tc corresponds to /c^ = 1/2. The 



above expressions are consistent with the following scal- 
ing form for the order parameter: 



(/)(z)=L-^/^X^(z/L,tLV^), 



(3.26) 



with /3 = u = 1/2. Note, however, that within the mean- 
field theory the magnitude of the scaling function is 
not universal, in that it is multiplied by the nonuniversal 
factor ^/2/u :^ 1.091. Finally, we stress that the choice 
of two parameterizations (see Eqs. ( 3.23) and_(3^25|) of 
the scaling functions in Eqs. (3.22) and (3.24), is just 



for convenience; it allows one to avoid using imaginary 
values of k and k. Indeed, one can transfer any of the set 
of equations into the other. For example, defining k as 



k 



where 



r / 2 



1 



(3.27) 



and taking into account the following properties of the 
elliptic functions [25l |26] 



K{k) = k' K{k), 



and 



dii{u] ik) 

sn(ii; ik sn(2xVTT^; k/^/l + k'^) ' 



(3.28) 



(3.29) 



3.22 



3.25 



one can easily check that the pair of equations 

(3.23) is equivalent to the pair of equations ( |3.24 ). ^ ^ 

In the this article — see Appendix |B] — we report the 
derivation of exact mean-field expressions for the be- 
havior of the local and of the total susceptibilities. We 
demonstrate that, when = 0, one has 



and 



X{xt) = L^X{xt), 



(3.30) 



(3.31) 



for the local and the total susceptibilities, respectively. 
The scaling functions of the total susceptibility X{xt) is 



X{xt) = 
where 



C2{xt)/K (k) ^ K (k) -2E jk) 



4.k'^k^K{k) 
k'^K(k)^(k^ -k'^)E(k)' 



(3.32) 



(3.33) 



Here E{k) is the complete elliptic integral of the second 
kind. The behavior of X{xt) is illustrated in Fig. [4] For 
the scaling function X^{z\xt) of the local susceptibility 
one has 



where 



¥"2 



1 - 2dn 



X, 



m,0 



C^k' 



(3.34) 



(3.35) 
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FIG. 4: Color online. The scaling function of the total sus- 
ceptibility X{xt). 



and 



U'2 



+/c cn I z / C; k 



k\l - 2k^')E ( am ( C; j ; j - ^^'X^,oC 

(3.36) 

i^'^ + (1 - 2^'^) dn (^i^C;k'^2 

with C e [-1/2,1/2] and 

X^,o = 2k'K{k) = 2K{k). (3.37) 
The behavior of the scaling function at T = Tc is 



0.015 



0.010 



0.005 



0.000 




gravity. The comparison is presented in Fig. [6j where 
Cr = {Sc^'^Kc)-^ =0.198. 
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FIG. 6: Color online. The scaling function of the total sus- 
ceptibility xL~^/^, calculated via the lattice model, compared 
versus the scaling function X{xt) of the same quantity as de- 
rived analytically within the continuum approach. 



IV. THE BEHAVIOR OF THE 
SUSCEPTIBILITY AT T = Tc 

Our analysis of the finite-size behavior of the system at 
T = Tc is summarized in Figs. [7| H [B EHl and[ll] Specif- 
ically, Fig. [7| presents the behavior of the normalized sus- 
ceptibility xL~^^^ in an experimentally realistic system 
in which both gravity and van der Waals interactions 
are taken into account, while Fig. [S] shows the behav- 
ior of the susceptibility in one idealized system in which 
only short-ranged interactions are taken into account. In 
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FIG. 5: Color online. Plot of the function X^{C\xt = 0). 
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shown on Fig. |5] 

As we will see in Appendix [Bj under proper rescaling 
of the abscissa and the vertical axes that follows from 
the mapping of the lattice onto the continuum model, 
there is perfect agreement between our lattice model re- 
sults for the total susceptibility and the analytically de- 
rived ones in the case of a short-ranges system without 



FIG. 7: Color online. The behavior of the finite-size suscep- 
tibility X at T = Tc as a function of Xp — /S.pL^^^ — 3.2 In L 
for films with thickness L = 1000, 2000, 4000 and 8000 layers. 
Both van der Waals forces and gravity effects are taken into 
account. 

both figures the behavior of the finite-size susceptibility 
is shown as a function of Xp = IS.pL^I'^ — clnL, where 
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c = 2.182 for the short-ranged systems and c = 3.2 for 
systems with van der Waals interaction present. Films 
with thickness L = 1000, 2000, 4000 and L = 8000 layers 
are considered, and 



1 ^ 



(4.1) 
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FIG. 8: Color online. The behavior of the finite-size suscepti- 
bility X at T = Tc as a function of Xp = ApL^^"" - 2.182 In L 
for films with thickness L = 1000, 2000,4000 and 8000 layers 
for systems with only short-ranged interaction present. 

The appearance of the InL corrections in the scal- 
ing variable X p IS a specific feature of the mean-field 
systems due to the degeneracy of the critical expo- 
nents /3 and u which become equal in mean-field ap- 
proximation. The numerical value of the constant c 
can be analytically predicted: at the bulk critical point 
K = Kc = {c2 ^ 2c^" )-^ - 0.168, h = one has a = 0, 
u = 1.681 and Eq. (3.17) possesses a solution, see Eq. 



(3.22), where the leading behavior near the boundary, 
when C ^ 0, is 



1.091/z. 



(4.2) 



Integrating over z and having in mind that the system 
is bounded by two substrate planes, one immediately ob- 
tains a InL contribution in Ap which is proportional to 
c = 2.182, i.e., exactly the same constant as given in 
Fig. m The presence of a van der Waals type interac- 
tion changes the constants of the model, e.g. the critical 
coupling Kc, as well as the effective constant in front 
of the second derivative of the order parameter profile. 
The last does not change the leading z-dependence of 
the order parameter profile but leads to a different con- 
stant c, which turns out to be c = 3.2 in our model. 
One observes that when both van der Waals forces and 
gravity effects are taken into account, see Fig|7| but L 
is not very large, the gravity effects are negligible - as 
for L = 1000,2000 and 4000; the corresponding curves 
are close to each other and one can speak about (some) 
data collapse for them. However, the curve for a system 



with L = 8000, for which the gravity effects are essen- 
tial, differs essentially from the others with the maximum 
of the susceptibility strongly suppressed. Furthermore, 
note that in the presence of van der Waals interactions 
and gravity the maximum of the curve for L = 8000 
shifts to higher values of Xp in comparison with systems 
with short-ranged interactions only and no gravity; see 
Fig. [Sj In the last case one observes a reasonable data 
collapse for all values of L considered. We stress that 
Ap contains contributions due to the role of the bound- 
aries which are not, in fact, critical. These are, e.g., the 
contributions which are due to the layers very near the 
boundaries, which layers are, independently on the value 
of Xp, always liquid like and almost fully occupied (i.e. 
densely packed) by the molecules of the fluid. Thus, in 
terms of Ap one expects larger corrections to scaling than 
when X is considered as a function of Ap (see Fig. [3|. 
The behavior of the finite-size susceptibility x at T = Tc 
as a function of Ap, as obtained within the mean-field 
like treatment of the model, is shown in Fig. [9j All the 
curves are for L = 1000 or L = 8000 layers thick film but 
with gravity and/or van der Waals interaction neglected 
or taken into account. One observes that both van der 
Waals interactions, as well as gravity, suppress the diver- 
gence of the susceptibility. Fig. [To] shows the same as 



Fig. [9] but now the behavior of the finite-size susceptibil- 
ity X is considered as a function of the scaling variable 
{/3Ap)L^^^ . One observes the different importance of 
the van der Waals substrate-fluid interactions and of the 
gravity within "thin" films, exemplified by L = 1000, and 
"thick" films, represented by the case L = 8000. While 
for L = 1000 the curves group together depending on 
the presence or absence of van der Waals substrate-fluid 
interactions, for L = 8000 they do this, at least with re- 
spect to the position of their maximum, depending on the 
presence of gravity in the system. Furthermore, one ob- 
serves that in "thick" films when ^ = the maximum of 
the susceptibility is at Ap < 0, contrary to the case with 
^ 7^ 0, when the maximum is at A/i > 0. In "thin" films 
the position of the maximum does not depend on g and 



always is at Ap < 0. Finally, Fig. 11 shows the density 
profiles of systems with L = 1000 and L = 8000 layers 
both for Xp = 0, as well as for the corresponding Xp for 
which the susceptibility has a maximum. More precisely, 
the part a) of this figure shows the density profile for a 
system with L = 1000 layers at T = Tc and p = pc. Note 
that Ap > everywhere, i.e., the equilibrium state of the 
film is liquid-like. For L = 1000 the effect of gravity is 
negligible, i.e., this profile is similar to the one with for 
L = 8000 when ^ = 0. Part b) of Fig. 11 shows the den- 
sity profile for a system with L = 1000 layers at T = Tc 
and Ap = Ap^ax for which the susceptibility reaches 
its maximum. Note that Apmax < 0. Next, part c) of 
Fig. [TT] presents the density profile for a system with 
L = 8000 layers at T = Tc and p = p^ Note that, due 
to the gravity, the profile is with Ap > near the walls, 
but Ap < everywhere in the middle of the system, i.e. 
in the middle of the system the equilibrium profile for 
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FIG. 9: Color online. The behavior of the finite-size susceptibility x at T = Tc as a function of Ap as obtained within the 
mean- field like treatment of the model. All the curves are for L = 1000 or L = 8000 layers thick film but with gravity and/or 
van der Waals interaction neglected or taken into account. One observes that both van der Waals interactions, as well as 
gravity, suppress the divergence of the susceptibility. 
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FIG. 10: Color online. The same as in Fig. [9] but now the behavior of the finite-size susceptibility x is considered as a function 
of the scaling variable {f3A/j,)L^^^. 



the finite film is ^a^-like. This has to be compared with 
case a) when Ap > for all < z < L. Finally, part 
d) of Fig. 11 presents the density profile for a system 
with L — 8000 layers at T = Tc and Ap = A/i^ax for 
which the susceptibility reaches its maximum. Note that, 
contrary to the case L = 1000, when the gravity is not 
important, A/imax > 0. 



V. THE BEHAVIOR OF THE SUSCEPTIBILITY 

FOR T>Tc AND T < Tc 

The behavior of the susceptibility away from the criti- 
cal point is summarized in Figs. [l2j [13] and [l4j More 
precisely. Fig. [12] presents the behavior of the finite- 
size susceptibility as a function of the scaling variable 

1000 and 
±10-^ When 
^ the system un- 



L 

t = 



ApL^I^ for a film with thicknesses L 
8000 for t = 0,t = ±10-^ and t -- 
-10~^ and at a negative value of x 



dergoes a second order phase transition for L = 1000, 
while for L = 8000 it undergoes first order phase tran- 
sitions for both t = — 10~^ and t = — 10~^. We further 
note that when L = 1000, the gravity effects being neg- 
ligible, the critical point Tc{L) of the finite system is at 
^IJ^c.L < 0, see Fig. [l] In general it is quite difficult to 
determine the location of this point with good precision 
Inspecting Fig. [12] one discovers, however, that for 
L = 1000 the curve with t = — 10~^ is very close to to 
the corresponding one that characterizes the behavior of 
the system at the true critical temperature of the finite 
system. In the context of the behavior displayed in this 



figure, it is useful to refer to Fig. [14] which illustrates the 
behavior of the excess normalized density Ap as a func- 
tion of the scaling variable = ApL^^^ for the same 
values of L and for the same fixed values of t as chosen 



in Fig. 12 



The behavior of the finite-size susceptibility as a func- 
tion of Ap for a film with thicknesses L = 1000 and for 
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FIG. 11: Color online, a) The density profile for a system with L = 1000 layers at T = Tc and fi = fic. b) The density profile 
for a system with L — 1000 layers at T = Tc and A/j = A/Xmax for which the susceptibility reaches its maximum, c) The 
density profile for a system with L = 8000 layers at T = Tc and fi = fic. d) The density profile for a system with L = 8000 
layers dX T — Tc and A/x = Afi^a^ for which the susceptibility reaches its maximum. In all the figures the bold (red) curve 
shows the profile of the local susceptibility normalized by for the system with size L. In order to have a better visibility of 
the order parameter profile simultaneously with the profile of the local susceptibility, the magnitude of the local susceptibility 
has been further reduced 8 times for L — 1000 and 16 times for L — 8000. 



t = 0, t = ±10-^ and t = ±10"^ is shown in Fig. [l3| Fi- 
nally, the behavior of Ap as a function of the scaling vari- 



A/iL^/^ for a film with thicknesses L = 1000 
8000 for^= 0,t = ±10-^ and t = ±10'^ 
Again, one observes that when 
^ the system undergoes a second 
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-10- 



able X ij 
and L 

is shown in Fig. 
L = 1000 and t = 
order phase transition at a negative value of while 
when L = 8000 it undergoes a first order transition for 
both t = — 10~^ and t = — 10~^ as is varied. We note 
that near Tc the position of the vapor-liquid coexistence 
line shifts (see Fig. [T]) when L increases, from a position 
characterized by A/i < to a position with Aja > 0. 
The same is also true for the position of the true criti- 
cal point of the finite system. For L = 8000 the critical 
point Tc{L = 8000) is at A/^^l > 0. The larger L the 
narrower the region in t and A/i in which one has round- 
ing of the second-order phase transition around the bulk 
critical point. Furthermore, we note that the larger L is 
stronger the gravity effects will be in the system; this in 
turn stimulates the phase separation within the system. 



These general remarks are intended to provide the 
reader with guidance regarding general trends of behav- 
ior as the system size L varies between 1000 and 8000. 
We also stress that the effects of van der Waals inter- 
actions, as well as those of gravity are not universal, in 
that they depend on the strength of two parameters: the 
value of the effective surface potential h^^g and the value 
of the effective gravity constant see Appendix [Aj Thus, 
detailed predictions for the finite-size behavior of the sus- 
ceptibility in a particular fluid bounded by specific sub- 
strates require that one perform numerical calculations 
following the general prescriptions presented in this arti- 
cle. 



VI. DISCUSSION AND CONCLUDING 
REMARKS 

We have studied the behavior of the finite-size suscep- 
tibility in fluid nonpolar films governed by van der Waals 
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FIG. 12: Color online. The behavior of the finite-size susceptibility as a function of the scaling variable = A/j.L^^^ for a 
film with thicknesses L = 1000 and L = 8000 for t = 0, t = ±10"^ and t = ±10"^ 
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FIG. 13: Color online. The behavior of the finite-size suscep- 
tibility as a function of Ap for a film with thickness L — 1000 
and for t = 0,t = ±10~^ and t = ±10~^ 



interactions and subjected to the influence of the gravi- 
tational field of the Earth. We focused on the situation 
in which the film is bounded by solid substrate plane 
boundaries which both strongly prefer the liquid phase 
of the fluid (i.e. we considered the so-called "plus-plus" 
boundary condition). 

Figure [15] summarizes the behavior of the susceptibility 
for such fllms as exemplifled with fllms with thicknesses 
L = 1000, 2000, 4000 and 8000 layers. The susceptibility 
is shown for temperature equal to that the bulk criti- 
cal temperature, T = Tc, plotted as a function of the 
scaling variable = A/j^L^^^ . The quantity A/i is the 
excess chemical potential where A/i > stabilizes the liq- 
uid phase of the fluid. Both the intrinsic van der Waals 
pair interactions between the molecules of the fluid and 
the van der Waals interaction between the fluid and the 
constituents of the two substrate plates are taken into 
account. For comparison, the behavior of the flnite-size 
susceptibility of a system with completely short-ranged 
interactions is also presented. We conclude that the be- 



havior of the susceptibility in realistic nonpolar fluid sys- 
tems, in which both van der Waals interactions and grav- 
ity are present, will always differ from the corresponding 
behavior of short-ranged systems, which constitute the 
standard theoretical model for such systems. As might 
have been expected, for realistic fllms with moderate fllm 
thickness L, the van der Waals interactions lead to no- 
ticeable differences in the behavior of the susceptibility 
from what one flnds for a short-ranged system. For large 
L, gravity gives rise to the the dominant effect. For the 
system studied here, consisting of He fllms bounded by 
Au surfaces it turns out that "moderate" thickness is a 
width of up to L = 4000 layers; the smaller the thickness 
the stronger the effect due to van der Waals interactions. 
It turns out that L = 8000 represents a "thick" fllm, for 
which the principal effects giving rise to the deviation 
from the short-ranged behavior of the susceptibility are 
due to the presence of a gravitational fleld in the region 
occupied by the fluid. At the bulk critical point T = Tc 
and A/i = for L = 8000 layers gravity gives rise to 
the appearance of a very large gas-like region that in- 
cludes the middle of the fllm which splits the two liquid- 
like regions near the solid substrate planes that bounds 
the fllm; see Fig. 11 This gas-like region occupies the 
greater part of both the upper and the lower half of the 
fluid system, but is larger near the upper bounding sur- 
face. This leads to the result that at T = Tc and A/i = 
and for L large enough, the average density of the flnite 
system is less than the critical density of the bulk system. 
Thus, in order to achieve coexistence between the fluid 
and the gas-like states of the system one needs to apply 
a positive excess chemical potential. This in turn leads 
to a new flnite-size coexistence line in thin fllms; see the 
topmost line in Fig. [l] which is different from the ones 
obtained on the basis of studies of fluid systems with no 
gravity present [T71 [El [El [2^ where such a shift of the 
coexistence line is always in the direction of negative A/i. 

Finally, we mention that in order to verify our lat- 
tice model approach we have in the case of a fully short- 
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FIG. 14: Color online. The behavior of as a function of the scaling variable = A/iL 
L = 1000 and L = 8000 for t = 0,t = ±10"^ and t = ±10"^ 



for a film with thicknesses 
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FIG. 15: Color online. Finite-size behavior of the susceptibil- 
ity for a film with L = 1000, 2000, 4000 and 8000 layers as a 
function of the scaling variable — A/iL^^^ in the case in 
which both the intrinsic van der Waals interactions as well as 
the presence of gravity are taken into account. For a compar- 
ison the behavior of a system with completely short-ranged 
interaction is also presented. 



ranged system derived analytically and compared with 
numerical calculations on a lattice model the behavior 
of the local and total susceptibilities. We find excellent 
agreement between the two models. The details are given 
in Appendix [B) 
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APPENDIX A: ESTIMATION OF THE GRAVITY 
CONSTANT IN TERMS OF THE PARAMETERS 
OF THE MODEL 



As we see from Eqs. (3.7) and (3.14), gravity intro- 



duces a field-like contribution into the grand canonical 
potential normalized per area. This contribution is re- 
fiected by the term 



1 ^ 



1 K 

2K~r 



Peg z(l){z) 



. (Al) 



Here we measure z from the bottom of the fiuid layer, 
which is positioned perpendicularly to the gravitational 
force. We work in units in which \(t){z)\ < 1 is a di- 
mensionless number, with ^(z) = 1 corresponding to 
one atom of ^He or ^He occupying a unit cell with vol- 
ume ttQ, where ao is the average distance between the 
helium atoms at the critical point. One can think of 
ao as being the lattice spacing of the lattice model. In 
[6j we estimated ao = 4.9 A for ^He and ao = 4.2 A 
for ^He; the basic material specific characteristics of the 
two isotopes of the helium needed for the current es- 
timations are summarized in Table [Tj Thus, the actual 
physical density, corresponding to (/)(z) = 1 is equal to 
p = uu/Oq^ where n is equal to 3 for ^He and 4 for ^He, 
and u is the atomic mass unit; lu = 1.6605 x 10~^^ kg. 
Next, with g = 9.81 m/s^ = 9.81J/(kg X m) and with 
z = /^aolO^^^m, where Iz is an integer that denotes 
the layer occupied by the corresponding particle is, one 
obtains 



Pcpgz = (2.2 X lO^^J-^) X 3 X 1.66054 x 10 



-27 



kg 



X 9.81 



0.527 X 10 



kg X m 

12 



X 4.9 X 10"^ V 



(A2) 
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for ^He, and 

Pcpgz = (1.4 X 10^^ J-^) X ^4 X 1.66054 x 10"^^^^ 

X ( 9.81 — I X X 4.2 X 10"^^ m) 
\ kg X m / ^ ^ 

= 0.383 X 10-^^-^ (A3) 
for ^He. In the above estimates we have taken into ac- 



TABLE I: Material-specific characteristics for ^He and ^He 
used to estimate parameters of the model investigated in the 
article. 





ao 


Tc 


/3c 


P 


Pc 


^He 


ao 4.9 A 


3.3 K 


2.2 X lO^^J"^ 


?>u/aQ 


0.041 g/cm^ 


^He 


ao 4.2 A 


5.2 K 


1.4 X lO^^J-^ 


Au/al 


0.069 g/cm^ 



effects. In [7] the conclusion has been made that we that 
gravity has a much more pronounced effect near liquid- 
gas critical points than for the ^He lambda transition. 
One can also compare the strength of the expected grav- 
ity effects in different substances near their respective 
critical points - see, e.g. [21] and Table I in [7 . There 
such a comparison is performed for ^He, SFe, Xe and 
CO2 in terms of the so-called "gravity scale height in a 
fluid" i^o, with Hq = Pc/{pcg), where Pc is the critical 
pressure of the fluid [21 . The smaller Hq the larger the 
gravity effect. It turns out that Hq is smallest for ^He, 
i.e. the gravity effect there is the strongest. Hq steadily 
increases for the sequence SFg, Xe and CO2, i.e. the 
strength of the gravity effects diminishes in these sub- 
stances in the sequence they are ordered in the current 
text. Our own estimates support this observation; we 
predict that the gravity effects in the critical behavior of 
^He will be stronger than in ^He, since the corresponding 
effective parameter reflecting the strength of the gravity 
in the model gs is larger than ^4; see above. 



count the fact that [6] 

f3c = 2.2 X 10^2 J"^ 

for ^He and 

Pc = 1.4 X 10^2 

for ^He (see Table [l|. We thus conclude that in our units 
one has to take the gravitational constant (times /3c) to 
be 

12 f_ 3t 



and 



{Peg) =g3 = 0.527 X 10"'^ for ''He 



{Peg) =g4 = 0.383 X 10"'^ for^He 



Note, howe ver, t hat in the equation for the order param- 
eter p rofile ( |3.1[ ), as well as in the excess grand potential 
( 3.14[ ), which are written in terms of (f) one always has 
a factor of 1/2 in front of gravitational constant. This 
should not be forgotten when performing a numerical 
evaluation of the gravity effect. 

Taking into account the fact that the critical part of the 
free energy near Tc behaves as with a = within 

the mean-field approximation, it is clear that gravita- 
tional effects can be felt in the thermodynamic behavior 
of the system when t'^ ~ 10~^^. This implies that one 
must explore relative temperature deviations of the order 
of t ~ 10~^ in order to be able to observe severe gravita- 
tional effects in the finite-size behavior of thermodynamic 
quantities. Obviously, the larger L, the stronger will be 
those effects. However, in such a case the observation of 
the finite-size properties of the studied quantities will be 
more challenging since they emerge when tL^/^ = 0(1). 
Thus, one needs to find a proper balance in the size of 
the system in order to observe both finite-size and gravity 



APPENDIX B: RESPONSE OF THE SYSTEM 
WITH SHORT RANGE 
INTERACTIONS^GINZBURG LANDAU 
APPROACH 

In the case in which all the interactions in the system 
are short-ranged and in the absence of gravity, Eq. (3.15) 



for the order parameter profile in a continuum system can 
be written in the standard form 



where 



1-^ 



and h ■. 



With the substitution 

4>iz) 

the above equation becomes 



m{z) 



am + 2m = h, 



where 



h = 



h. 



The equations for the layer response function 



d4>{z) 



dh 



dm{z) 



h=0 



dh 



(Bl) 
h 

cfK' 
(B2) 

(B3) 
(B4) 
(B5) 

(B6) 



h=Q 
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then reads 



+ (d + 6m^)xi = 1- 



Introducing the scahng variables 



Xt = aL 



(B7) 



(B8) 



one can rewrite the equations (B4) and (B7) for the order 



parameter m{z) and of the local susceptibility Xi{^) i^^o 
equations for the scaling functions 



and 



Xm{z\xt, Xh) = L m{z\xt, Xh) 



X^{z\xt,Xh) = L xi{A^t,Xh) 



of these quantities. One obtains 
(fx, 



XtXm + Xh, 



and 



^ + (a:* + 6X^)X^ = l, 



(B9) 



(BIO) 



(Bll) 



(B12) 



respectively. Because boundary conditions are identi- 
cal at both surfaces of the film system, the solutions 
of the above equations have to satisfy (j)'{L/2) = and 
x'{L/2) = or, equivalently, X'^iC = 1/2) = and 
X!^{( = 1/2) = 0, i.e., the middle of the system will be 
an inflection point for all physical quantities which de- 
pend on the distance from the boundaries. Note that 
despi te th e fact that there is no explicit dependence of 
Eq. (B12) on x^, the scaling function, X^, do es de pend 
on Xh because X^, which is a solution of Eq. ( Bll|), de - 
pends on that parameter and Xm enters into Eq. (B12). 



In the remainder of this appendix we will be interested 
in the behavior of X^ on z and Xt for Xh = 0. Then, the 



solution X m{z\x t) of Eq. (Bll) is known. From Eqs. 
(|3.22|) and (|3.24|) one has 



a) when Xt > — tt^ 

where < C ^ 1^ /c G M is to be determined 
form Eq. ([3^. 



b) when Xt < — tt^ 



2K{k) 



sn[2K{k)C;~ky 



(B14) 



where < ^ ^ 1^ and /c G M is to be determined 
formEq. ([3^. 




FIG. 16: Color online. Plot of the profile function Xm{C\^t) 
- see Eq. \B17\. 



In order to utilize the symmetry of the problem it is 
helpful to move the coordinate frame so that the origin 
of the system is at the midpoint of the film. Taking into 
account that [24 



and that 



we obtain 



d-n[u^K{k)]k\ _ k' 
sn[u -\- K{k)]k] cn(ii; /c) ' 

cii{iu]k') = l/cn(ii;/c). 



Xra{C\xt) = Xrafi Cn\i2K {k)C, k'] 

where C e [-1/2,1/2] and 

Xm,o = 2k'K{k) = 2K{k). 
Since cn(0; k) = 1, one has 

Xmfi = Xm{0\xt)- 



(B15) 

(B16) 
(B17) 
(B18) 
(B19) 



A typical behavior of X^{C,\xt) is shown on Fig. [Tg] 

Given the scaling function Xm all that one has to do 
to determine X^{C,\xt) is to solve Eq. (B12) with the 
boundary conditions X^(±l/2) = and dX^{() / d( = 
for C = 0. 

According to the general theory of differential equa- 
tions of second order 

X^{C\xt) = ciMCl^t) + C2MCM + c,^^{C\xt), (B20) 

where ci, C2 and q are constants, ipi and ip2 are linearly 
independent solutions of the homogeneous equation 



d^X^ 
dC 



(x, + 6X^)X^ = 0, 



(B21) 



and ij ji is a particular solution of the inhomogeneous Eq. 
(B12). It is easy to check that 



MC\xt) = ^X^{C\xt) = X^{C\xt) (B22) 
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FIG. 17: Color online. Plot of the function - see Eq. 

(B23l. 



is a solution of Eq. (B21). Explicitly, one has 



^i{C\xt) = -^^sn [i2K{k)C; k^] dn [i2K{k)a k'] 



-I — -^sn 



dn 



.(B23) 



A typical behavior of the solution of the homogeneous 
equation ipi{(\xt) is shown on Fig. 17 Because of the 
symmetry of the problem, 



Since, see Eq. (|B23|, 



(B24) 
(B25) 



one concludes that ci = 0; see Eq. (B20). 




FIG. 18: Color online. Plot of the function ip2{C\^t) - see Eq. 
(B29l. 



Following Abel [271 [28], we construct V^2(Ckt) 

dC 



(B26) 



Obviously, '02(Cl^t) is an even function with respect to 
(. Again, because of the symmetry of this system, one 



has dilj2{C\xt)/d( = at (" = 0. Next, it is easy to check 
that the Wronskian of the solutions of tpi and ?/^2 

w.Mc/-^-Mc/-^ (B27) 

is equal to one, i.e., the two solutions i/ji and 1^2 are lin- 
early independent. Explicitly, performing the integration 
for ip2^ one obtains 



fc'(l - 2fc'2)E ( am ( (^■k'];k']- ik^X^^oC 



k 



+A.'cn(^^C;A.') 



(B28) 

X fc'2 + (i_2fc'2)dn(^i^C;fc')2 |. 
A typical behavior of the solution of the homogeneous 



equation V^2(Ckt) is shown on Fig. 18 
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FIG. 19: Color online. Plot of the function ij^i^Qxt) - see Eq. 
(B3TI. 



Again following Abel [271 HF ^ one finds the particu- 
lar solution of the inhomogeneous equation expressed in 
terms of and ^1)2 via 

m\xt) = Mc\xt) I M<:\xt)d<:-M<:\xt) J Mc\xt)dc. 

(B29) 



Taking into account Eq. {B22), one can show that 

X„,{C\xt) 



V'i(Ckt) = -V': 



(B30) 



Using this result and that fact that t/ji is a solution of 
the homogeneous equation (B21), one can show that ipi 
is indeed a solution of Eq. ( |B12 ). Furthermore, since 
the constant on the right-hand side of Eq. ( B12| ) is equal 
to one, we have q = 1. One can explicitly determine 
i^iiCl^t)- Performing the integration in Eq. (B30), one 
obtains 



U/2 



1 - 2dn 



X 



m,0 



(B31) 
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with dtlji{(\xt)/d( = at C = 0. A typical behavior 
of the particular solution of the inhomogeneous equa- 
tion tlji{(\xt) is shown on Fig. 19 Since q = 1 
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FIG. 20: Color online. Plot of the function xi{0- 



from X^{l/2\xt = 0) we determine that 
-lim^^i/2[^i(Ckt)/V^2(Cl^t)]- Explicitly, after 



and ci 

C2{Xt) = 

series of manipulations one can show that 



(B32) 



The behavior of the scaling function T = Tc 

obtained from the continuum approach described above 
is shown on Fig. |5]in the main text. One can also calcu- 
late the local susceptibility profile from the lattice model 
approach. The result for L = 1000 is shown on Fig. 



20 One can easily check that the curves in both fig- 



ures completely overlap each other if L~'^xu calculated 
via the lattice model, is multiplied by = 0.212. The 
last is very close to what one theoretically predicts from 



the mapping of the lattice model onto the continuum 
one. Inspection shows that = {?>C2^Kc)~^ = 0.198. 
The above demonstrates that the lattice model correctly 
reproduces the universal scaling function in the case in 
which that function can be calculated exactly and, thus, 
is trustworthy for predicting properties of the local and 
total susceptibility in cases when no analytical results are 
available. Such an instance, considered in the main text, 
is one of a system in which both gravity and van der 
Waals interactions are present. 

In addition to the scaling function of the local suscep- 
tibility one can also derive in an analytic closed form the 
scaling function X of the total susceptibility. Starting 
from 



X{xt) 



/•1/2 

7 ■ 

Jo 



X^{C\xt)dC, 



one derives 



X{xt) 



c2{xt)/K {k) ^ K {k) -2E (k) 



(B33) 



(B34) 



The behavior of this function is illustrated in Fig. [4] in 
the main text. 

In order to compare with the corresponding result from 
the lattice calculations one must take into account the 
fact that, due to the difference in definitions of the field 
variable, X must be rescaled by 1/c^. Additionally, tak- 
ing in to a ccoun t the fact that {l — K/Kc)L^^^ = CrXt (see 
Eqs. ( |B2| ) and (B8)) one obtains the result shown in Fig. 
|6] in the main text. We are thus able to conclude that 
the lattice model produces a result in a p erfect agreement 
with the analytic solution in Eq. (B34) for the behavior 
of the total susceptibility. 
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